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We study the Renyi entropy of the one-dimensional XYZ spin-1/2 chain in the entirety of its 
phase diagram. The model has several quantum critical lines corresponding to rotated XX Z chains 
in their paramagnetic phase, and four tricritical points where these phases join. Two of these points 
are described by a conformal field theory and close to them the entropy scales as the logarithm of its 
mass gap. The other two points are not conformal and the entropy has a peculiar singular behavior in 
their neighbors, characteristic of an essential singularity. At these nonconformal points the model 
undergoes a discontinuous transition, with a level crossing in the ground state and a quadratic 
excitation spectrum. We propose the entropy as an efficient tool to determine the discontinuous or 
continuous nature of a phase transition also in more complicated models. 
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In the past several years there has been a constantly 
increasing interest in quantifying and studying the en- 
tanglement of virtually every physical system p~r[3j . This 
interest is not surprising, as entanglement provides valu- 
able insights from many different perspectives. 

As it is often measured as entanglement entropy, that 
is, the Von Neumann entropy of the reduced density ma- 
trix of a system S = — Trplnp, its origin lies in the con- 
text of quantum information theory 4, 5 . Essentially, it 
quantifies the "quantumness" of a state and therefore it 
provides a measure of its suitability for efficient quantum 
algorithms, in the ongoing quest for quantum computing. 

In the physics of strongly interacting systems, entan- 
glement has been welcomed as a new interesting correla- 
tion function, different in nature compared to the tradi- 
tional ones, due to its nonlocal structure. In particular, 
it has provided a new challenge for the integrable model 
community on one side 6 8 and has led to new, more ef- 
ficient approaches for numerical simulations (tensor net- 
work states) on the other [91-11 1 j . For statistical physics, 
the entanglement entropy has also been proposed as a nu- 
merically efficient way to characterize a system, due to 
its diverging behavior across a phase transition [T^l H3] ■ 

If we concentrate on the bipartite entanglement, that 
is, the entanglement entropy between two complementary 
regions, a lot is understood about its qualitative behav- 
ior. In general, it satisfies the so-called area law [2], that 
is, it is proportional to the area of the boundary dividing 
the two regions. This is naively understood considering 
that the entanglement is carried by correlations, that, 
in general, decay exponentially with the distance with a 
rate given by the correlation length £. If the volume of 
the two regions is much bigger than £, entanglement is 
localized at the boundary. Hence the area law. This be- 



havior is modified when correlations decay more slowly, 
like close to phase transitions, where £ — > 00. 

The one-dimensional (ID) case is most interesting. 
As the boundary between regions consists of individual 
points, in a massive phase the entropy saturates to a finite 
value for sufficiently large systems. For critical phases 
described by a conformal theory, correlations decay as 
power laws and the entropy diverges logarithmically with 
the size of the region [T3j. Close to conformal points, the 
behavior is still governed by a Conformal Field Theory 
(CFT) and the divergence is logarithmical in £ [15] . 

In this letter, we argue that close to a nonconformal 
point of phase transition the entropy shows a dramati- 
cally different behavior, characteristic of an essential sin- 
gularity. Thus, we propose the entropy as an efficient 
(analytical or numerical) tool to distinguish, for exam- 
ple, a first-order ferromagnetic Quantum Phase Transi- 
tion (QPT) in a spin chain from higher-orders ones. This 
kind of ability can be very helpful in the study of a vari- 
ety of other models with nonconformal critical (but still 
gapless) points, such as those described by spin-1 Hamil- 
tonians |16j or in fermionic systems |17j and many more. 

An essential singularity for the entropy was already 
observed at the bicritical point of the XY model in [18] . 
where it was noticed that the conformal prediction fails 
because low-energy excitations have a quadratic disper- 
sion relation. However, the XY chain can be considered 
a toy model; peculiar behaviors like these can be simply 
a mathematical feature, that does not survive in more re- 
alistic systems. Here we will consider a fully interacting 
model that cannot be mapped into free systems, like the 
XY chain and those of ,22J, that is, the XYZ spin-1/2 
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chain, denned by the Hamiltonian 

HxYZ = - ^2 { J x a nK+l + J V a n a n+l + J * a n a n+l) > 
n 

(i) 

where <7„ (a = x, y, z) are the Pauli matrices acting on 
the site n, the sum ranges over all sites n of the chain 
and the constants J x , J y and J z take into account the 
degree of anisotropy of the model. Without any loss of 
generality, we can rescale the energy and set J x = 1, 
J y = T, and J z = A. This model is the most general 
spin-1/2 ID chain with nearest neighbor interactions and 
its rich phase diagram has both conformal and noncon- 
formal critical points. Moreover, unlike |18j . this is a 
truly interacting, although still integrable, model. 

Following the standard procedure, we start from the 
ground-state wavefunction | ) of and divide the sys- 
tem into two parts, which we take as the semi-infinite 
left and right chains, thus dividing the Hilbert space as 
H = Hr^Hl- We introduce the reduced density matrix 
by tracing out one of the two half-chains, say T-Ll'- 



p = Tr HL | 0><0 



(2) 



As an entanglement estimator we will consider the Renyi 
entropy [19] 



1 — a 



lnTr/3" 



(3) 



which reduces to the Von Neumann entropy for a — > 1 
(a being just an additional parameter at our disposal). 

In the thermodynamic limit, the reduced density ma- 
trix of the XYZ model is known exactly, due to its con- 
nection with the Corner Transfer Matrices (CTM) of the 
zero-field 8- vertex model [20] . Using this connection, the 
entropy of the XYZ chain was calculated in [2"T] . where 
it was shown that p can be written in the form 



1 

Z 
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(4) 



where x = exp[— 7rA//(fc)], Z = + an d k and 

A are elliptic parametrizations of the coupling constants 
that, in the principal regime of the 8-vertex model, can 
be written as 



l + fcsn 2 (iA) 
1 - k sn 2 (iA) ' 



cn(iA) dn(iA) 
~l-/ssn 2 (iA) 



(5) 



with 0<fc<l, 0<A< I(k'). I(k) is the complete 
elliptic integral of the first kind and k' = yl — k 2 is the 
conjugated elliptic parameter. 

In [22] , the entanglement was studied for systems akin 
to free particles, using the fact that p is of the form Q 
with eigenvalues Xj = e -2 -? 6 or Xj = e -1 - 2 ^ 1 - 16 , for the 
ordered or disordered phases, with e characteristic of the 
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FIG. 1: Phase Diagram of the XYZ model in the (A, V) plane. 
The solid lines -A = ±1, |r| < 1; T = ±1, |A| < 1 and 
A = ±r, |A| > 1- correspond to the critical phase of a ro- 
tated XXZ chain. Out of the four "tricritical" points, Ci t % 
are conformal and E\ 2 are not. 



model. From the CTM of [20] it is known that all inte- 
grable, local spin-1/2 chains have a reduced density ma- 
trix of the form Q. Thus, the degeneracy of the eigen- 
values of p is universal for integrable models (and given 
by a partitioning problem |23| ) and the entanglement is 
completely characterized by e. In the case of the XY Z 



model, 6 — it 27^y • 

Further defining q = e~ 



a 



a 



the Renyi entropy is 

l + T^X>( 1 + ^ Q ) ■ 

1 — a ^— ' v ' 



which we can rewrite in terms of elliptic 9 functions 
as: 



S a {e) = 



a 



6(1 



In 



In 



t(0,g)g 3 (0,g) 
2 2 (0,q a ) 



(6) 



(7) 



ln(2) 
3 



6(1 -a) 9 3 {0,q a )9 4 (0,q«) 

This form, although completely equivalent to ([6]), is most 
suitable for studying the behavior of the entropy in the 
phase diagram of the model, thanks to the vast literature 
on elliptic functions. 

In fig. [T] we plot a graph of the phase diagram of the 
XYZ chain. The model is symmetric under reflections 
along the diagonals in the (A,T) plane. The system is 
gapped in the whole plane, except for six critical half- 
lines/segments: A = ±1, |r| < 1; T = ±1, |A| < 1 and 
A = ±r, |A| > 1. All of these lines represent param- 
agnetic XXZ chains, with the anisotropy along different 
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FIG. 2: Contour plot of the Von Neumann entropy of the 
XYZ model in the (A,F) plane. Regions of similar colors 
have similar entropy values and the lines where colors change 
are lines of constant entropy. The brighter is the color, the 
bigger is the entropy. 



directions. They are thus described by a c = 1 (ultravi- 
olet) sine-Gordon theory with /3 dependent on A or T. 
These critical lines meet three by three at four points, 
that we will denote as "tricritical" , since they separate 
three regions where the system is gapped. These three 
regions can be obtained one from the others by exchang- 
ing the role of the three spin components a x , a y , and 
o z . Two of these points -C\^ = (1, — 1), ( — 1, 1)— are 
conformal points with f3 2 — 8tt: while the other two - 
Ei_ 2 = (1, 1), (— 1, — 1)- correspond to j3 = and are 
nonconformal. The former points correspond to an an- 
tifcrromagnetic Heisenberg chain at the BKT transition, 
while the latter correspond to a Heisenberg ferromagnet. 
Thus, at Ei 2 the system undergoes a transition in which 
the ground state passes from a disordered state to a fully 
aligned one. Exactly at the transition, the ground state 
is highly degenerate while the low-energy excitations are 
magnons with a quadratic dispersion relation. 

The different behavior of the entropy close to these 
points is clear from fig. [2j where we show a contour plot 
of the Von Neumann entropy as a function of T and A. 
In any neighbor of the entropy is diverging, since 
it is controlled by a conformal theory. Close to E\ 2 the 
entropy goes from to infinity depending on the direction 
of approach to the nonconformal point. Moreover, we see 
that all curves of constant entropy pass through either 
Ei t 2- This is the same behavior observed around the 
bicritical point of the XY model in [TB] . 



The nonconformal point of the latter model does not 
show a highly degenerate ground state, however, the den- 
sity of states of low-energy excitations always diverges 
like n(e) ~ for a quadratic dispersion e(k) ~ k 2 . 

Thus, in both models the ground state close to the multi- 
critical point can be constructed out of a thermodynami- 
cally large number of configurations. This argument pro- 
vides a qualitative picture of how this rich variety in the 
entanglement may arise close to nonconformal points. 

For the XYZ model we can make these statements 
more quantitative by expanding the parameters around 
the tricritical points. Close to G2 we take 

r = 1 - <5cos0 , A = -l-<Ssin0, < < | , 

yielding 

k = (Jtim <j) + 1 - Vtan + O (5) , (8) 

a = sj ^ (v/wTT - Vw) + o (<5 3 / 2 ) . 

Thus, k is not defined at the tricritical point, since its 
value as 5 — > depends on the direction <ft of approach 
(for (f) = 0, k = 1 and for (f> = 7r/2, k = 0), while A — > as 
5 — > 0, irrespective of 0. Hence, e ~ \/S and in any neigh- 
borhood of the conformal point the entropy is diverging. 
In particular, on the direction = 0, since 1(1) — > 00, 
the entropy is divergent for any S (as expected, since this 
is a critical line). 

Expanding instead around E\ as 

r = -l + <Scos0, A = -l-<5sin0, < <j> < ^ , 

we find 

k = ( Vtan (j) + 1 - Vtan ^ + O (8) , (9) 

sn(iA) = u - 2 , (Vtan0+ 1 - J tan A + 0(<5 1/2 ) . 
V cos <p \ J 

Thus, k is exactly as in the previous case, while A has 
a very singular behavior. A regular expansion of A in 
powers of <5 is not readily available. However, one can 
see that A ~ I(k') close to E\ and thus e ~ I(k')/I(k), 
which varies from for <f> = to infinity for cf> = tt/2. 

As we have already mentioned, the behavior of the 
entropy is completely determined by e, thus, the Renyi 
entropy has an essential singularity at E\^ and can take 
any positive real value arbitrarily close to it. This un- 
usual behavior can be understood qualitatively consider- 
ing that this tricritical point is a highly symmetric one 
with a highly degenerate ground state. Thus, as the pa- 
rameters are changed to pass through £1,2, the ground 
state undergoes a level crossing and the entanglement 
can change in a discontinuous way. This effect is most 
dramatic along one of the XX Z lines, where the entropy 
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goes suddenly from diverging in the paramagnetic phase 
to vanishing in the ferromagnetic one. 

From the examples of the XYZ chain and of the XY 
model |18j . we can infer that points of discontinuous 
phase transitions show a characteristic singular behavior 
of the entanglement entropy, which can be easily detected 
numerically and might be used as an efficient tool to 
distinguish, for example, between first- and higher-order 
phase transitions in more complicated models [TBJ IT7] . 
It is important to remark that S a and the correlation 
length show very different behaviors close to i?i,2, since, 
as it is expected for any phase transition, £ diverges as 
one approaches the critical lines or points, without show- 
ing the essential singularity characteristic of the entropy 
[2"5] , Thus, close to the discontinuous points the en- 
tropy formula derived in |14j is no longer valid: indeed it 
applies in the vicinity of all conformal points of the criti- 
cal lines (including the BKT points Ci^), while it fails if 
we are close to nonconformal points, such as £4.2, where 
the spectrum of excitations has a quadratic dispersion 
relation for small momenta. 

The behavior of the Renyi entropy close to a conformal 
point has been recently investigated in |15j . where the 
Poisson summation formula was applied to ^ to obtain 
its asymptotic behavior. Eq. ^ is equally suitable for 
such an expansion, using the different representations of 
the Jacobi theta functions In [T5], the Ising model 
and the XXZ chain were considered. Using the XYZ 
and the elliptic properties, we can easily extend their 
analysis. For instance, we can expand ^ in powers of 
t = e _7r / e for e -> close to G\ 2: 



S a (e) 



t z 



(10) 



t+- + 0(t 3 ) 



J-2/lc 

t 1 ^ + -^ + (i 3/Q ) 



(l-a) 
1 

"(T^o) 

Because, close to these two tricritical points we have 
r 1 =s 4e- 2 / 2e = At 1 / 2 , 



we can write 



S a 
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12 



1 



Hi) 



a 

(l^a) 
1 



i- 2 
16 
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o(r 



-2/a 



-Mi) 



-4/a 



o i 



-6/a 



We recognize the familiar leading logarithmic term, a 
nonuniversal constant, and a series of corrections, in 
agreement with the results of [T5]. For large a the lead- 
ing correction is given by with x = 2. This im- 
plies that the operator responsible for this correction has 



conformal dimensions (A, A) = (1,1) and is therefore a 
quadratic combination of the SU(2)i Kac-Moody cur- 
rents J, J. 

Conclusions We have studied the bipartite Renyi en- 
tropy of the ID XY Z spin chain in its phase diagram, 
using exact analytical expressions derived from the in- 
tegrability of the model. The entropy diverges on the 
critical lines: close to conformal points the divergence 
is logarithmical in the correlation length with power-law 
corrections. At the nonconformal points the entropy has 
an essential singularity. We argued that this may be a 
characteristic feature of discontinuous phase transition 
points that could allow to easily numerically discrimi- 
nate between first- and higher-order phase transitions. It 
would be interesting to test this behavior in other models. 

Expression ^ for the Renyi entropy is general for inte- 
grable models of spin 1/2 and possesses very interesting 
elliptic and modular properties, as a function of k, A and 
a. Different modular properties, in real space in- 

stead of parameter space, emerge close to the conformal 
points and have been used in |15j to understand the lead- 
ing corrections to the entanglement entropy. It would be 
interesting to understand whether there is a link between 
these two. 
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